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The Astrophysical Journal; in press 
ABSTRACT 

We have used the 28' x 28' Hubble Space Telescope image mosaic from the GEMS (Galaxy Evolution from 
Morphology and SEDs) survey in conjunction with the COMBO-17 photometric redshift survey to constrain the 
incidence of major mergers between spheroid-dominated galaxies with little cold gas (dry mergers) since z = 0.7. 
A set of yV-body merger simulations was used to explore the morphological signatures of such interactions: they 
are recognizable either as < 5 kpc separation close pairs or because of broad, low surface brightness tidal features 
and asymmetries. Data with the depth and resolution of GEMS are sensitive to dry mergers between galaxies with 
My < -20.5 for z < 0.7; dry mergers at higher redshifts are not easily recovered in single-orbit HST imaging. 
Six dry mergers (12 galaxies) with luminosity ratios between 1:1 and 4:1 were found from a sample of 379 red 
early-type galaxies with My < -20.5 and 0.1 < z < 0.7. The simulations suggest that the morphological signatures 
of dry merging are visible for ~ 150 Myr and we use this timescale to convert the observed merger incidence into 
a rate. On this basis we find that present day spheroidal galaxies with My < -20.5 on average have undergone 
between 0.5 and 2 major dry mergers since z ~ 0.7. 

We have compared this result with the predictions of a Cold Dark Matter based semi-analytic galaxy formation 
model. The model reproduces the observed declining major merger fraction of bright galaxies and the space 
density of luminous early-type galaxies reasonably well. The predicted dry merger fraction is consistent with our 
observational result. Hence, hierarchical models predict and observations now show that major dry mergers are 
an important driver of the evolution of massive early-type galaxies in recent epochs. 

Subject headings: galaxies: general — galaxies: interactions — galaxies: elliptical and lenticular — galaxies: 
evolution — galaxies: structure 



1. introduction 

Numerical simulations have long predicted that early-type 
galaxies, with spheroid-dominated stellar light profiles, are 
a natural outcome of major galaxy mergers (e.g., Toomre & 
Toomre 1972; Barnes & Hernquist 1996; Naab & Burkert 2003). 
Yet, both the timing and nature of the violent assembly of early- 
type galaxies remain frustratingly unclear. Look-back studies 
have recently become large enough to demonstrate a steady 
growth in the total stellar mass in the red-sequence galaxy pop- 
ulation since z ~ 1 (e.g., Chen et al. 2003; Bell et al. 2004a; 
Faber et al. 2005). The majority of these red-sequence galax- 
ies are morphologically early-type out to at least z ~ 0.7 (Bell 
et al. 2004b); indeed, recent works have confirmed a growth 
in the total stellar mass in morphologically early-type galaxies 
from z ~ 1 to the present day (Cross et al. 2004; Conselice et 
al. 2005). This growth is dominated by a growing number of 
~ 0.5 -2L* early-type galaxies (e.g., Bell et al. 2004a; Drory et 



al. 2004; Faber et al. 2005) 1 . Interestingly, relatively few blue 
galaxies bright enough to be the star-forming progenitor of a 
massive non-star-forming early-type galaxy are observed (Bell 
et al. 2004a; Faber et al. 2005). Thus, it has been suggested that 
mergers between non-star-forming early-type galaxies may oc- 
cur (so-called 'dry mergers'), and build up the massive early- 
type galaxy population (Bell et al. 2004a; Faber et al. 2005). 
Such mergers are observed at least in clustered environments 
(van Dokkum et al. 1999; Tran et al. 2005). Yet, it is still un- 
clear how important dry mergers are in driving the evolution of 
luminous, massive early-type galaxies, averaged over all cos- 
mic environments. 

Dry mergers may also play an important role in explain- 
ing another interesting phenomenon: the strong differences be- 
tween the observed properties of high-luminosity (with My < 
-21) and low-luminosity early-type galaxies. Luminous (i.e., 
massive) early-type galaxies possess boxy isophotes, cores in 

1 Although it is possible that the number of very massive galaxies, with L > 
5L*, does not change by a large amount since z ~ 1-5 (e.g., Saracco et al. 2005). 
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their light distributions, and are supported primarily by ran- 
dom stellar motions. In contrast, faint early-type galaxies have 
diskier isophotes, derive more support from organized rota- 
tional motions, and have power-law (or Sersic) light distribu- 
tions all the way into their central-most regions (e.g., Bender 
1988; Kormendy & Bender 1996; Lauer et al. 1995; Faber et al. 
1997; Graham 2004). These structural differences may be the 
results of galaxy merging between different types of progeni- 
tors (Khochfar & Burkert 2003). 

There is growing evidence that the faint early-type galaxy 
population has grown substantially at recent times through 
galaxy merging and/or disk fading. There is a deficit of faint 
red-sequence galaxies in z ~ 1 galaxy clusters (Kodama et al. 
2004; de Lucia et al. 2004a), and those that are present appear to 
have rather young stellar populations (van der Wei et al. 2004). 
Furthermore, faint early-type galaxies at the present day can 
have important contributions from younger stars (Kuntschner 
2000; Trager et al. 2000; Thomas et al. 2005). Simulations of 
merging disks easily yield remnants with properties consistent 
with the typically disky, partially rotationally-supported low- 
luminosity elliptical galaxies (e.g., Toomre & Toomre 1972; 
Naab, Burkert, & Hernquist 1999; Bendo & Barnes 2000; 
Naab & Burkert 2003), especially for unequal-mass interac- 
tions. Taken together with observations of a declining merger 
rate, implying between 0.2 and 1 merger per L* galaxy (e.g., 
Le Fevre et al. 2000; Patton et al. 2002; Conselice et al. 2003, 
although see Lin et al. 2005), it is likely that mergers between 
gas-rich galaxies have played an important role in building up 
low-luminosity early-type galaxies. 

In contrast, it is increasingly clear that mergers of gas-rich 
disks cannot be the main formation route of more massive 
early-type galaxies. Not only are there relatively few blue 
star-forming galaxies at z < 1 luminous enough to fade into 
present-day massive non-star-forming galaxies, but the boxy 
isophotes and cores in the light distributions of massive early- 
type galaxies are challenging to reproduce in mergers of disk 
galaxies (Naab & Burkert 2001; Barnes 2002; Naab & Burkert 
2003). Such properties can, however, be naturally reproduced 
by mergers of elliptical galaxies (Naab, Khochfar, & Burkert 
2006). This hypothesis is given some observational support by 
van Dokkum et al. (1999) and Tran et al. (2005), who found 
that ~ 15% of bright galaxies in the outskirts of the cluster 
MS 1054-03 were in bound pairs and were likely to merge on 
< 1 Gyr timescales; most of these are pairs of early-type galax- 
ies. Semi-analytic work by Khochfar & Burkert (2003, 2005) 
supports the formation of massive early-type galaxies ending 
up in clusters by this route. 

A direct measure of the frequency of dry mergers over the 
last half of cosmic time will further constrain the importance of 
this process in driving the evolution of the luminous early-type 
galaxy population. Furthermore, because the majority of mas- 
sive galaxies are red early-type galaxies at all epochs since z ~ 1 
(Bell et al. 2004a), the dry merger rate gives important insight 
into the merging history of the most massive galaxies at recent 
times. Finally, gas-free mergers are relatively straightforward 
to model, and their luminosity and color evolution is easily pre- 
dicted, significantly simplifying the estimation of merger rates. 

The goal of this paper is to use these data to constrain the 
incidence of dry mergers in the GEMS dataset. To realize 
this goal, a number of ingredients must be in place. Firstly, 
a scheme must be devised by which one can reliably and re- 
producibly identify mergers of gas-poor systems. Secondly, the 



timescale over which signatures of gas-poor merging are iden- 
tifiable must be characterized. Finally, these criteria must be 
applied to the data; coupled with the timescale, a merger fre- 
quency can be estimated. 

In this paper, we present our first attempt to address this issue 
using high-resolution HST data from GEMS (Galaxy Evolu- 
tion from Morphology and SEDs; Rix et al. 2004), photometric 
redshifts from COMBO- 17 (Classifying Objects by Medium- 
Band Observations in 17 Filters; Wolf et al. 2004), and state- 
of-the-art A^-body and semi-analytic simulations. We briefly 
describe the data in §2. We outline the selection philosophy 
in §3. In order to devise criteria by which one can select can- 
didate dry mergers and understand over which timescale these 
selection criteria are sensitive, we explore the characteristics 
of a suite of A^-body simulation mergers between spheroid- 
dominated galaxies (§3.1). We then apply these selection cri- 
teria to the data in §3.2. We use these results to derive a dry 
merger fraction and incidence in §4. In §5 we compare our 
results with limits on dry merger rates from other observations, 
compare with published merger rates, and compare the observa- 
tions with the predictions of a Cold Dark Matter (CDM) based 
semi-analytic model. In the Appendix, we explore the impor- 
tance of high-speed fly-by (unbound) interactions, and argue 
that the vast majority of tidally-distorted interacting spheroids 
should be in the lowest-density environments capable of host- 
ing spheroid-dominated galaxies, and are therefore likely to 
merge. Throughout, we assume f2 m = 0.3, fi m + f2A = 1, and 
H = 70kms- 1 Mpc- 1 . 

2. THE DATA 

COMBO- 17 has imaged the ~ 30' x 30' extended Chandra 
Deep Field South using 5 broad and 12 medium passbands com- 
plete to apparent /?-band magnitude limits of niR ~ 23.5. Using 
these 17-passband photometric data in conjunction with galaxy, 
star, and AGN template spectra, classifications and redshifts are 
assigned for ~ 99% of the objects with sufficient flux. The 
typical galaxy redshift accuracy is &/(l +z) ~ 0.02 (Wolf et 
al. 2004), allowing construction of rest-frame colors and abso- 
lute magnitudes accurate to ~ 0.1 mag. In order to probe inter- 
actions between gas-free, non-star-forming galaxies, we select 
only galaxies on the red sequence with photometric redshifts 
0.1 < z < 0.7 for further study, following Bell et al. (2004a) 
and Mcintosh et al. (2005). 

We use F850LP imaging from the GEMS survey to provide 
0"07 resolution data for our sample of red-selected galaxies. 
Using the Advanced Camera for Surveys (Ford et al. 2003) on 
the HST, a ~ 28' x 28' area of the Extended Chandra Deep Field 
was surveyed to a depth allowing galaxy detection to a limiting 
surface brightness of ^f850lp,ab ~ 24 mag arcsec" 2 (Rix et al. 
2004); in practice early-type galaxies have relatively high sur- 
face brightnesses, and are limited only by the magnitude limit 
of COMBO- 17 redshift classification. 

3. THE SELECTION OF DRY MERGERS 

For the sake of clarity, we present our selection criteria for 
candidate dry mergers here. Our philosophy is to catch the 
merging galaxies during the early stage of merging, where the 
two galaxy nuclei are distinct entities allowing relatively accu- 
rate characterization of the progenitors, including, importantly, 
an estimate of the mass ratio of the merger. 

Therefore, we select for further study very close pairs of 
red-sequence bulge-dominated galaxies with i) projected sep- 
aration < 20kpc, ii) a V-band absolute magnitude difference of 
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Table 1 

Properties of the dry merger candidates 



Galaxy name 


RA 


Dec 


z 


M v 


U-V 


Classification 




hh mm ss 


dd mm ss 










Galaxies thought not to be merging 


GEMS 033139.25m275303.0 


03 31 39.25 


-27 53 03.0 


0.20 


-21.5 


1.1 


00000 


GEMS 033139.03m275300.1 


03 31 39.03 


-27 53 00.1 


0.20 


-21.4 


1.3 


00000 


GEMS 033145.09m275520.6 


03 31 45.09 


-27 55 20.6 


0.49 


-20.9 


0.9 


03003 


GEMS 033145.09m275517.4 


03 31 45.09 


-27 55 17.4 


0.57 


-20.0 


1.0 


03003 


GEMS 033155.67m273840.0 


03 31 55.67 


-27 38 40.0 


0.519 


-20.8 


1.2 


03300 


GEMS 033155.83m273839.4 


03 31 55.83 


-27 38 39.4 


0.526 


-20.1 


1.2 


03300 


GEMS 033307.89m275 150.4 


03 33 07.89 


-27 51 50.4 


0.518 


-21.6 


1.3 


00003 


GEMS 033307.75m275151.2 


03 33 07.75 


-27 51 51.2 


0.44 


-20.3 


1.0 


00003 


Contentious cases 


GEMS 0331 17.77m2758 12.1 


03 11 17.77 


-27 58 12.1 


0.58 


-21.0 


1.0 


30303 


GEMS 033117.82m275814.1 


03 11 17.82 


-27 58 14.1 


0.55 


-20.3 


1.0 


30303 


GEMS 033313.59m275735.3 


03 33 13.59 


-27 57 35.3 


0.66 


-23 


1.25 


20310 


GEMS 033313.51m275737.3 


03 33 13.51 


-27 57 37.3 


0.65 


-22.5 


1.1 


20310 


GEMS 033123.29m274544.9A 


03 31 23.29 


-27 45 44.9 


0.55 


-20.8 


0.95 


22101 


GEMS 033123.29m274544.9B 






0.55 


-19.3 


0.95 


22101 


Mergers of galaxies with significant bulges and gas 


GEMS 033244.06m275 143.3 


03 32 44.06 


-27 51 43.3 


0.26 


-21.7 


1.3 


33333 


GEMS 033244.27m275 141.1 


03 32 44.27 


-27 51 41.1 


0.273 


-20.5 


0.95 


33333 


GEMS 033215.97m274943.1 


03 32 15.97 


-27 49 43.1 


0.59 1 


-22.3 


1.5 


33033 


GEMS 033216.16m274941.6 


03 32 16.16 


-27 49 41.6 


0.67 


-22.1 


0.9 


33033 


GEMS 033151.21m274659.3 


03 31 51.21 


-27 46 59.3 


0.67 2 


-22.7 


1.0 


23333 


GEMS 033151.37m274700.3 


03 31 51.37 


-27 47 00.3 


0.64 


-21.9 


1.2 


23333 


Galaxies classified as dry mergers by the majority of classifiers 


GEMS 033136.86m280128.0A 


03 31 36.86 


-28 01 28.0 


0.15 


-22.1 


1.5 


23222 


GEMS 033136.86m280128.0B 






0.15 


-20.6 


1.5 


23222 


GEMS 033157.38m274239.0A 


03 31 57.38 


-27 42 39.0 


0.62 


-22 


1.3 


22222 


GEMS 033157.38m274239.0B 






0.62 


-22 


1.3 


22222 


GEMS 033208. 85m275851.8 3 


03 32 08.85 


-27 58 51.8 


0.54 


-20.4 


1.0 


23321 


GEMS 033208.76m275851.2 3 


03 32 08.76 


-27 58 51.2 


0.54 


-20.5 


1.0 


23321 


GEMS 033210.77m274234.4A 4 


03 32 10.77 


-27 42 34.4 


0.419 


-21.9 


1.2 


22222 


GEMS 033210.77m274234.4B 4 






0.419 


-21.9 


1.2 


22222 


GEMS 033239.48m275301.6 


03 32 39.48 


-27 53 01.6 


0.65 


-22.5 


1.25 


22221 


GEMS 033239.47m275300.5 


03 32 39.47 


-27 53 00.5 


0.65 


-21.8 


1.2 


22221 


GEMS 033300.99m275545.2A 5 


03 33 00.99 


-27 55 45.2 


0.620 


-21.7 


1.0 


23321 


GEMS 033300.99m275545.2B 






0.620 


-21.7 


1.0 


23321 



Note. — 1 The redshift difference between GEMS 033215.97m274943.1 and GEMS 033216.16m274941.6 is somewhat larger than expected for COMBO-17's 
typical redshift error; either one of the redshifts is substantially in error or our visual classification is inappropriate. 

2 The redshift difference between GEMS 033151.37m274700.3 and GEMS 033151.21m274659.3 is consistent with COMBO-17's redshift error (the galaxies' visual 
morphologies show that they are interacting, and are therefore at the same redshift). 

3 This merger candidate is one COMBO-17 object but was successfully split by the GEMS source extraction algorithm. The coordinates are taken from GEMS, the 
redshift and color are adopted directly from COMBO-17, and the absolute magnitude of the COMBO-17 object is split into both galaxies using the observed F850LP 
magnitude ratio from GEMS. 

4 GEMS 033210.77m274234.4A and B have been recently discussed in a Keck laser guide star adaptive optics study of mergers detected by Chandra (Melbourne et 
al. 2005, their XID-536; with an X-ray luminosity of nearly 10 42 ergs/sec). Their luminosities and redshifts are consistent with our own. They fit model SEDs to the 
observed 0.4-2.2^tm broad band fluxes, concluding that both components are dominated by old stellar populations at the epoch of observation. 

5 GEMS 033300.99m275545.2A and B have a COMBO-17 redshift of 0.45; the VVDS redshift is adopted in this paper. Absolute magnitudes and rest-frame colors 
have been re-calculated using the correct VVDS redshift. 

GEMS galaxy names are simply a combination of their RA and Dec (where 'm' denotes negative declinations), and A and B refer to the two sub-components when 
the system was not resolved into two separate galaxies by COMBO-17. Rest-frame My values are approximate, split using our rough estimate of the luminosity ratio 
of the merger. Redshifts are accurate to <5z ~ 0.02, magnitudes to ~ 0.2 mag, and colors to ~ 0.15 mag. Redshifts quoted to 3 significant figures are spectroscopic 
redshifts from the VVDS (Le Fevre et al. 2004). See the text for details of the classification scheme. 
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FIG. 1 . — A simulation sequence showing the final coalescence of a gas-free elliptical-elliptical 1:1 merger. This particular simulation had an initially parabolic 
orbit with a pericenter of 21 kpc; the face-on projection of the orbital plane is shown in this example. The merger remnant has a final M v = -22.5 and is placed at 
z = 0.5 in real background noise from GEMS with the right size and apparent magnitude. In the first frame, the orbital paths of each galaxy are shown (ticks are 
every 26 Myr), where a time of corresponds to first pass, 2.21 Gyr to second pass, and final coalescence happens at 2.48 Gyr. The total elapsed time between the 
first and last frame is 340 Myr, starting 2. 15 Gyr after first pass as the galaxies are coming towards second pass (the third frame) and final coalescence (the last frame 
of the sequence). The images are 50 kpc on a side. This interaction is classifiable as a merger for roughly 200 Myr via either faint and broad tidal features or very 
close double nuclei. 



1 .5 mag or less (corresponding to a 4: 1 luminosity ratio differ- 
ence, our threshold for a major merger), Hi) photometric red- 
shift difference |Az| < 0.1. 

However, not all projected close pairs will be merging; some 
may be projections of galaxies at significantly different red- 
shifts, and some may be unbound fly-by interactions which will 
not merge. In order to address these problems, it was necessary 
to add a second layer of selection. Prolonged experimentation 
with automated merger diagnostics did not yield satisfactory re- 
sults, and were unable to robustly flag the existence of broad, 
low surface brightness sheets of debris which are the hallmark 
of a gas-poor major merger. Therefore, we have used visual 
classification to screen the automatically-selected subsample of 
close pairs for signs of dry merging. This second level of vi- 
sual selection is very stringent: in demanding low-level tidal 
features or very close separations, many merging systems are 
thrown out of the sample in phases where the tidal features 
are unobservably weak (simulations suggest this should be the 
case around ~ 1/2 of the time). While this reduces the sample 
size, it provides some measure of security against projections 
and fly-bys, and results in a smaller but cleaner sample of dry 
merger candidates (spectroscopy of our dry merger candidates 
is being sought and will be discussed in future works). 

In order to calibrate out and understand the subjectivity of 
visual classifications, we have adopted a dual strategy: we have 
visually classified a suite of simulated dry mergers and fly- 
bys(selected to have projected separation < 20 kpc) in order to 
understand better the signatures of a dry merger and to char- 
acterize the timescale over which the signs of dry merging are 
recognizable (§3.1), and we use 5 independent sets of classifi- 



cations to quantify the degree of reproducibility of dry merger 
classifications (§3.2). 

3.1. Expected Morphological Signatures of Dry Mergers 

In order to explore the expected appearance of morphological 
indicators of dry mergers, and the timescales over which these 
indicators are visible with single-orbit HST imaging, we make 
use of the simulations of Naab, Khochfar, & Burkert (2006). 
The progenitors are early-type galaxies, and are themselves the 
remnants of gas-free disk-disk mergers from Naab & Burkert 
(2003) 2 . There are 400000 particles (160000 stars and 240000 
dark) per progenitor. A number of simulations with different 
pericenters and orbital parameters (initially marginally bound, 
parabolic, or a hyperbolic fly-by) were carried out for two mass 
ratios: 1:1 and 3:1. A variety of different viewing angles were 
chosen between face-on and edge-on (with respect to the or- 
bital plane). The simulation output was recorded every 26 Myr. 
The simulations were used to construct artificial GEMS im- 
ages; these were created by smoothing to GEMS resolution, re- 
binning to the appropriate angular diameter distance, and scal- 
ing to reproduce a given total luminosity of the remnant. The 
images were then added to real sky background from GEMS. 

From examination of these simulated images (e.g., Fig. 1), 
some important points become clear. The morphological sig- 

2 Gas-rich major merger remnants are somewhat diskier than remnants from 
gas-free major disk mergers (Barnes 2002, Naab & Burkert, in prep.). The 
adoption of a gas-rich major merger remnant as the spheroid merger progenitor 
in what follows does not significantly affect the outcome, while avoiding ex- 
tra uncertainties associated with the detailed choice of star formation rate and 
feedback prescriptions. 
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FIG. 2. — Six simulation snapshots deemed to be non-interacting by the majority of the classifiers. We show the simulation and the GALFIT residual side-by-side 
for each postage stamp. Images are 5()kpc on a side. The residuals in the centers of the galaxies are from small mismatches between the PSF used to broaden the 
simulations and the GEMS PSF. 




FIG. 3. — Six simulation snapshots deemed to be interacting by the majority of the classifiers. The format is similar to Fig. 2. In all but the top left and bottom 
right postage stamps, the merging system was determined by SExtractor to be one object and was therefore fit as one object by GALFIT. There is further discussion 
of this issue in the text. 



natures of interactions between early-type galaxies are often 
very weak; tidal tails are broad with low surface brightness ow- 
ing to the high velocity dispersion of the progenitor. Further- 
more, there are nearly order-of-magnitude variations in overall 
merger timescales (i.e., the time from first pass to final coales- 
cence), which depend sensitively on parameters such as mass 
ratio, pericenter distance and orbit. Importantly, however, the 
timescale over which one can recognize tidal features in dry 
mergers is largely independent of these parameters, and is com- 
parable with the galaxies' internal dynamical times. Tidal fea- 
tures of sufficient strength to be detected in our images occur 
only when the galaxies are very close together — typically be- 
tween the last close pass prior to coalescence and coalescence 
itself 3 . Orientation plays a role: interactions viewed nearly par- 
allel to the orbital plane are visible as very close pairs for 2- 
3 times longer than the same interactions viewed from < 60° 
from face-on. 

3 In simulations with small pericenter distances tidal features are visible for 
a short time after first pass, but become rapidly too weak to observe. 



In order to estimate timescales over which GEMS-depth data 
will be sensitive to dry mergers, we explored the visibility of a 
wide range of early-type galaxy interactions at a number of red- 
shifts, and assuming a range in final merger remnant luminos- 
ity. Preliminary explorations showed that luminous dry merg- 
ers at z > 0.7 are visible only as a very close pair (timescales 

< 50Myr, except for rare projections along the orbital plane); 
the data lack sufficient depth to robustly detect tidal features. 
Accordingly, we limit this study to z < 0.7 galaxies and simu- 
lations. 

Simulation snapshots were selected randomly from a vari- 
ety of interaction sequences; the only selection criterion ap- 
plied was that the nuclei of the galaxies should be separated by 

< 20kpc (this is an analogous selection to that applied to the 
data). A total of 174 snapshots from 15 simulations with dif- 
ferent mass ratios (1:1 and 3:1), orbital parameters (marginally 
bound, parabolic, and a hyperbolic fly-by) and viewing angles 
(viewing angles from 0° from the normal to the orbital plane, 
30°, 60°, 80°, and the edge-on 90° orientation) met this sepa- 
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FIG. 4. — The four candidates which were not deemed to be interacting by the majority of the classifiers. We show the image and the GALFIT residual side-by-side 
for each sample system. Images are 50 kpc on a side. 




FIG. 5. — Three contentious candidates, arranged in order of decreasing star formation as judged from the classifications. GEMS 033123. 29m274544.9 was 
deemed by SExtractor to be one object, and was fit by GALFIT by a single Sersic profile. 



ration criterion 4 . The random selection was weighted to ensure 
that the different viewing angles were selected approximately 
uniformly, in terms of solid angle. Given that a simulation snap- 
shot is created every 26 Myr, this implies that the simulations 
have two distinct nuclei separated by < 20 kpc for on average 
300 Myr. The snapshots were scaled to a range of luminosities 
My < -20.5 and redshifts 0.1 < z < 0.7 5 , and a subsample of 
25 snapshots were SExtracted 6 , GALFITted and classified by 
the five classifiers (EFB, TN, DHM, CW and SK) in exactly 
the same way as the data. The classification was blind: simula- 
tion snapshots were intermixed with real dry merger candidates 
from GEMS. 

Systems were deemed to be not merging if 2 classifiers or 
fewer classified the system as a merger, possibly merging if 3 
classifiers argued it was a merger, and merging if 4 or 5 clas- 
sifiers decided it was a merging system. According to this 
scheme, 12/25 snapshots were deemed not to be merging and 
13/25 were declared merging (there were no cases in which 
only 3 classifiers declared a system to be merging). Examples 
of each class are shown in Figs. 2 and 3. It is apparent that 
snapshots with very weak/nonexistent tidal features were typi- 

4 Only one of the fly-by snapshots had separation < 20 kpc, reflecting the 
rapidity of fly-bys. 

5 The redshifts were drawn from a uniform deviate in terms of volume, not 
redshift, and are therefore weighted towards higher redshift in exactly the same 
way as the observational dataset (the average redshift of the simulations was 
0.53, whereas that of the observational dataset is 0.5). 

6 Some very close pairs were deemed to be one object by SExtractor and 
were therefore fit with a single Sersic profile. 



cally classified as non-interacting, while snapshots with either 
i) < 5 kpc separation and/or ii) tidal distortions (tails or asym- 
metries) were classified as interacting. Accordingly, we assign 
a timescale of ~ 150 ± 50 Myr over which dry mergers are both 
automatically selected (with separations < 20 kpc; timescales 
~ 300 Myr) and visually classified as a merger (13/25 of the 
snapshots, plus associated counting uncertainties). 

3.2. Selection of Dry Merger Candidates 

The dry merger simulations indicated that dry merger iden- 
tification is a well-posed problem with quantifiable selection 
effects for galaxies with z < 0.7 and My < -20.5: at larger red- 
shift and/or fainter limits tidal features are no longer easily vis- 
ible in single-orbit depth HST/ACS imaging. Accordingly, we 
adopt as the parent sample the sample of red sequence galaxies 
with U-V > 1 - 0.3 lz-0.08(M v + 20.77), 0.1 < z < 0.7, and 
My < -20.5 — a sample of 468 galaxies. Galaxy morpholo- 
gies were assigned by eye for these galaxies: see e.g., Mcin- 
tosh et al. (2005), Wolf et al. (2005), or Bell et al. (2004b) for 
a description of the morphological typing criteria and for some 
examples from each type. A total of 379/468 galaxies are clas- 
sified to have early morphological types, E/S0/Sa 7 . 

Adopting the 468-galaxy red sequence 0.1 < z < 0.7 and 
M v < -20.5 sample as the parent sample, we have selected all 
close pairs satisfying the following criteria: ;) projected sepa- 
ration < 20 kpc, a V-band absolute magnitude difference of 

7 Systems with small amounts of ongoing star formation in a weak disk com- 
ponent were allowed into the sample, as long as they were on the red sequence 
and their light distributions were dominated by a smooth bulge component 
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FIG. 6. — Three interactions with clear signs of dust and/or star formation, indicating some gas content. In all cases, both interaction partners have colors consistent 
with the red sequence (despite oftentimes small amounts of star formation in these systems), and substantial bulge components. 



1 iE.B6rnZB0 1 2B.0 






it | 






GEW5r033MD.ggm275£45.2 

. ; 1 


■' ■ 
. ... m%L. ■ ■ 



FIG. 7. — The six candidate spheroid-dominated merging systems in GEMS; the last 2 systems were judged by 2 out of the 5 classifiers to have hints of 
ongoing star formation. Images are 50kpc on a side. The spiral galaxy in the lower left-hand side of the GEMS 033136.86m280128.0 A and B postage stamp is a 
background projection. GEMS 033136.86m280128.0, GEMS 033157.38m274239.0 GEMS 033210.77m274234.4, and GEMS 033300.99m275545.2 were deemed 
by SExtractor to be single objects, and were therefore fit by GALFIT using single Sersic profiles. 



1 .5 mag or less (corresponding to a 4: 1 luminosity ratio differ- 
ence, our threshold for a major merger), Hi) photometric red- 
shift difference |Az| < 0.1. Ten pairs were found with these 
properties. 

However, from the visual classifications there were a num- 
ber of merging systems which were classified by COMBO- 17 
as one object. It is important to augment the sample of merger 
candidates with such blended mergers. Every one of the 468- 
galaxy parent sample has been fit using the GALFIT galaxy 
fitting code Peng et al. (2002) using a single Sersic profile 8 . 
A sample of 34 galaxies was automatically selected which had 
strong residuals from the Sersic profile fit (within one half-light 
radius, the sum of the absolute values of the residuals exceeded 
15% of the total flux within the half-light radius). The vast ma- 
jority of these 34 systems were bulge/disk systems or had strong 
dust features and lanes: only 6 systems had evidence for mul- 
tiple nuclei. These 6 systems were added to the sample of 10 
COMBO- 17-selected close galaxy pairs, forming the final sam- 
ple of 16 dry merger candidates. Postage stamps and GALFIT 
residuals of all 16 candidates can be found in Figs. 4-7. 



8 The Sersic (1968) model has a profile with surface brightness £ oc e~ r " , 
where r is the radius and n is an index denoting how concentrated the light 
distribution of a given galaxy is: n = 1 corresponds to an exponential light 
distribution, and n = 4 corresponds to the well-known de Vaucouleurs profile. 



4. RESULTS 

The 16 dry merger candidates were classified by the 5 classi- 
fiers (EFB, TN, DHM, CW and SK), with 4 different classifica- 
tion categories: (0) no evidence for an interaction, (1) possible 
merger, (2) dry merger, and (3) gas-rich merger. The results 
from all classifiers are given in Table 1 . On the whole, agree- 
ment between the classifiers was good, and the galaxies sepa- 
rated reasonably cleanly into 4 subsamples: galaxies where the 
majority of classifiers felt that there was no evidence for merg- 
ing (Fig. 4), systems where classification was uncertain (Fig. 5), 
gas-rich mergers (Fig. 6), and dry mergers (Fig. 7). A firm clas- 
sification as a merging system required that 4 or 5 classifiers felt 
that a given system was a definite merger. A total of 6 systems 
were identified as dry mergers, three were identified as gas-rich 
mergers (these systems all had prominent bulge components 
and red colors), and 2 of the possible mergers could qualify as a 
dry merger as both components were bulge-dominated (GEMS 
033123.29m274544.9 and GEMS 033313.59m275735.3). Fig. 
8 shows the color-magnitude diagram of all 0.1 < z < 0.7 
galaxies from COMBO- 17 (points), all red-sequence early-type 
galaxies (gray circles), and the merger candidates (black sym- 
bols, where symbol type denotes classification). 

On this basis, it is possible to estimate the incidence of ma- 
jor mergers between luminous early-type galaxies from z ~ 0.7 
to the present day. Among the 379 red, early-type galaxies in 
GEMS with M v < -20.5 and 0. 1 < z < 0.7, we find a total of 12 
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galaxies (6 pairs) that are strong dry major merger candidates. 
Classifications of dry merger simulation snapshots suggest that 
we will be able to recognize a dry merger for roughly an internal 
dynamical time, or approximately 150 Myr. Assuming a con- 
stant dry merger fraction over the 6.3 Gyr interval < z < 0.7, 
an average, luminous M v < -20.5 early-type galaxy will un- 
dergo (6.3/0.15) x (12/379) - 1.3^ dry major merger (with 
mass ratios between 1 : 1 and 4: 1). Yet, because there are only 6 
merging pairs, we cannot strongly constrain the merger fraction 
evolution. The evolution of the merger fraction has been argued 
to evolve as rapidly as (1 +z) 3 (e.g., Le Fevre et al. 2000). If the 
dry merger fraction were to evolve similarly, an average lumi- 
nous My < -20.5 early-type galaxy will undergo <~ 0.9±q J ma- 
jor dry merger between z = 0.7 and the present day. Obviously, 
our limited number statistics are currently a dominant source 
of uncertainty; nonetheless, our results indicate that dry merg- 
ing plays an important role in shaping the properties of lumi- 
nous early-type galaxies and building up the mass contributed 
by red-sequence galaxies since z ~ 1 . 

5. DISCUSSION 
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FIG. 8. — The distribution of the candidate spheroid mergers in the color- 
magnitude diagram. Points denote all galaxies with 0.1 < z, < 0.7; grey circles 
early-type galaxies on the red sequence. Black symbols show possible mergers 
(crosses), gas-rich mergers (diamonds), and dry mergers (filled circles). In the 
cases where COMBO- 17 was unable to resolve the merger from the ground 
and assigned a single redshift, we split the flux between the two galaxies fol- 
lowing Table 1 and assign the same color; small random offsets are added to 
aid visibility. The arrow shows the approximate position of L* for red se- 
quence galaxies at z ~ 0.5, and the dotted line denotes the locus of galaxies of 
constant stellar mass ~3x 10 10 Mq, assuming a Kroupa (2001) stellar initial 
mass function. 



5.1. Are our observations consistent with other observational 

constraints? 

Is our suggestion that nearly every luminous early-type galaxy 
in the nearby universe has experienced a major gas-free merger 
since z ~ 1 consistent with existing observational constraints? 
Dissipationless (=gas-free) merging preserves the fundamental 
plane (e.g., Nipoti, Londrillo, & Ciotti 2003; Gonzalez-Garcfa 
& van Albada 2003; Boylan-Kolchin, Ma & Quataert 2005); 
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FIG. 9. — A comparison of the dry merger fraction with previous observa- 
tions of merger fractions, and the predictions of a semi-analytic model. Gray 
data points show the merger fraction of galaxies with Mg < -20 from Patton 
et al. (2002, open circles), Le Fevre et al. (2000, open diamonds), Conselice 
et al. (2003, solid circles), Lin et al. (2005, triangles), and Bundy et al. (2004, 
crosses, and the upper limit at z = 0.63). The model prediction of the major 
merger rate for Mg < —20 galaxies is shown as a solid gray line for the red- 
shift range 0.1 < z. < 1.1. The black data point shows the dry merger fraction 
from this work. The black naked error bar shows the model prediction for the 
fraction of Mv < -20.5 red early-type galaxies that have experienced a major 
merger in the last 150 Myr averaged over 0. 1 < z < 0.7 with associated number 
uncertainties. 



therefore the thickness of the fundamental plane is not a strong 
constraint on the prevalence of dry merging. In contrast, sim- 
ulations and models have demonstrated that dry merging will 
lead to scatter in the stellar mass-size relation (equivalently, 
the Kormendy relation; Nipoti, Londrillo, & Ciotti 2003) and 
the color-magnitude relation (CMR; Bower, Kodama, & Ter- 
levich 1998). Observations of the evolution and scatter of these 
scaling relations limit the amount of dry merging to around a 
factor of two increase in mass since z ~ 1 (Mcintosh et al. 2005; 
Bower, Kodama, & Terlevich 1998) 9 . A similar limit is derived 
from the stellar mass evacuated from the cores of luminous el- 
liptical galaxies by their merging black holes in the absence of 
gas (Graham 2004). 

5.2. Comparison with previous measurements of merger 
fraction 

It is interesting to compare our results with previous obser- 
vational estimates of the overall galaxy merger fraction. In 
Fig. 9, we show the dry merger fraction (black solid point) 
of My < -20.5 red early-type galaxies, along with the frac- 
tion of M B < -20 galaxies in close pairs with < 20 kpc sepa- 
ration taken from Le Fevre et al. (2000, open diamonds) and 
Patton et al. (2002, open circles). The results for Lin et al. 
(2005) are shown as triangles (separations between 10 and 
30 kpc; M B < -20). Bundy et al. (2004) argued that many close 
pairs with nearly equal optical magnitudes were minor merg- 
ers where the satellite galaxies had boosted optical flux from 
tidally-induced bursts of star formation; they used near-infrared 

9 The CMR limit applies since the epoch when the CMR was imprinted on 
the early-type galaxy population. 
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photometry to estimate the major merger fraction (crosses, and 
the upper limit at z = 0.63). The solid circles denote the merger 
fraction of Mb < -20 galaxies, inferred using morphological 
disturbance as an indicator of past merging activity (Conselice 
et al. 2003). 

It is clear that the dry merger fraction is lower than the gen- 
eral merger fraction. A number of factors may contribute to 
this reduction. Close pairs and morphologically-selected gas- 
rich mergers are expected to be visible for a timescale r ~ 
0.5-1 Gyr (Le Fevre et al. 2000; Patton et al. 2002; Conselice 
et al. 2003), whereas dry mergers are visible for fd ry ~ 150 Myr. 
Episodes of star formation in gas-rich mergers will increase the 
brightness of the merging galaxies, increasing their prominence 
relative to the star formation-free dry mergers. Finally, gas-rich 
galaxies tend to inhabit lower-density environments than gas- 
poor galaxies; this difference in environment will also affect 
the likelihood of a major merger. 

5.3. Theoretical Expectations in a Hierarchical Universe 

Do hierarchical models predict that dry mergers are an im- 
portant mode of transformation and mass growth for luminous 
elliptical galaxies? To address this question, we make use of an 
updated version of the Somerville et al. semi-analytic galaxy 
formation model (see Somerville & Primack 1999; Somerville, 
Primack, & Faber 2001, for a description of the basic ingre- 
dients), based in the standard ACDM paradigm of hierarchical 
structure formation. The models are based on Monte Carlo re- 
alizations of dark matter halo merger histories, constructed us- 
ing the method described in Somerville & Kolatt (1999), and 
include fairly standard recipes treating gas cooling, star for- 
mation, and supernova feedback. The star formation rate is 
given by m* = m C oid/(Tj?fd y nX where m co id is the cold gas mass 
in the galaxy, fdyn is the dynamical time of the disk, and r° is 
a free parameter, which is adjusted to obtain agreement with 
observed gas fractions in spiral galaxies at the present day. We 
convolve the resulting star formation histories with the multi- 
metallicity stellar SED models of Bruzual & Chariot (2003), 
and include the effects of dust extinction using a simple rela- 
tionship between face-on B-band optical depth and star forma- 
tion rate (jb oc SFR 13 ), based on the observational results of 
Wang & Heckman (1996, we assume (3 = 0.5). Merging of 
galaxies (sub-structure) within virialized dark matter halos is 
tracked by computing the time it takes for a satellite galaxy to 
lose all of its orbital angular momentum via dynamical friction 
against the background of the dark matter halo, using the stan- 
dard Chandrasekhar approximation. Mergers with mass ratio 
1:4 or greater are assumed to result in a burst of star forma- 
tion, and all pre-existing 'disk' stars are transferred to a 'bulge' 
component. 

Earlier versions of this model, like many other models in 
the literature, had difficulty producing enough luminous red 
galaxies, particularly at high redshift (see, e.g., Somerville et 
al. 2004), and therefore could not produce a useful prediction 
of the merger rate for these systems. Somerville et al. (in 
prep.) find that introducing an ad hoc 'quenching' of star for- 
mation when the mass of the bulge component grows larger 
than Mbuige.crit — 2 x 10 10 M Q leads to much improved agree- 
ment with the observed < z < 1 color distributions as a func- 
tion of absolute magnitude, including the observed bimodal- 
ity of the color distributions at z ~ 0-2. A full description of 
this model will be presented in this forthcoming work: at this 
stage, we argue that as the model approximately reproduces the 



global luminosity functions and observed number of red, early 
type galaxies at z ~ 0-1, it is at least a useful starting point to 
examine the overall merger rate and the dry merger rate. 

In Fig. 9, we show the predicted fraction of Mb < -20 galax- 
ies that have experienced a major merger within the past 1 Gyr 
(the grey line; major mergers are defined as having mass ra- 
tios between 1:1 and 4:1). These predictions are derived from 
a mock light cone from 0.1 < z < 1.1 covering a sky area of 
1 sq. deg. (corresponding to about four times the volume of 
GEMS). Bearing in mind on one hand the difficulties of mea- 
suring galaxy merger rate from close pairs, and on the other 
hand uncertainties in physical prescriptions for dynamical fric- 
tion, tidal stripping of galaxy halos, and tidally-induced star 
formation, the agreement between the model and observations 
is rather good. There is a hint that the model merger fraction de- 
creases rather less rapidly than the observations; interestingly, 
this behavior is also seen in the semi-analytic models of Benson 
et al. (2002), as shown in Fig. 13 of Conselice et al. (2003). 

The black naked error bar shows the dry merger fraction pre- 
dicted by this model. We identify dry mergers in the model 
as galaxies with magnitudes My < -20.5, U — V colors on the 
red sequence, and which have had a major merger in the past 
150 Myr. The mock catalog has roughly four times the volume 
of GEMS; there are 13 merging systems from 982 early-type 
M v < -20.5 red sequence galaxies in this mock catalog in the 
interval 0.1 < z < 0.7. The inferred merger fraction is in ex- 
cellent quantitative agreement with our observational determi- 
nation (6 systems from 379 galaxies). This result lends fur- 
ther theoretical support to hydrodynamical studies of elliptical 
galaxy formation (Dormnguez-Tenreiro, Saiz, & Serna 2004) 
and semi-analytic studies of galaxy clusters (Khochfar & Burk- 
ert 2003), both of which argue that dry mergers are an important 
process in driving the buildup of massive early-type galaxies at 
recent times. 

6. CONCLUSIONS 

We have used the GEMS and COMBO- 17 surveys in con- 
junction with TY-body and semi-analytic galaxy formation sim- 
ulations to explore the frequency of gas-free major mergers be- 
tween spheroid-dominated galaxies (dry mergers) since z = 0.7. 
We focused on this aspect of galaxy merging both because iden- 
tifying gas-free mergers is less likely to be complicated by the 
effects of merger-induced star formation and the accompany- 
ing dust extinction, and because it can be modeled relatively 
straightforwardly (for similar reasons). The morphological sig- 
natures of such interactions were calibrated using mock GEMS- 
like images based on A^-body simulations, and include < 5 kpc 
separation of close nuclei and/or broad tidal tails and asymme- 
tries. The simulations showed that such features are visible only 
for galaxies with M v < -20.5 in single-orbit HST F850LP data 
at z < 0.7; higher-redshift galaxies show unobservably weak 
tidal features in the GEMS data. A total of 809 red-sequence 
galaxies were visually-inspected for morphological signatures 
of spheroid merging; 6 systems (12 galaxies) with luminosity 
ratios between 1 : 1 and 4: 1 were found out of a total of 379 red 
sequence early-type galaxies with My < -20.5. The simula- 
tions suggest that morphological signatures of spheroid merg- 
ing are visible for ~ 150 Myr; we therefore argue that an aver- 
age My < -20.5 early-type galaxy has experienced between 0.5 
and 2 major spheroid mergers since z ~ 0.7. This merger fre- 
quency is consistent with limits of < 1 major gas-free merger 
at recent times, as derived from the color-magnitude and stel- 
lar mass-size relations, and from core mass deficits. We have 
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compared this result with an updated semi-analytic model. This 
model reproduces the evolution of the number of luminous red 
early-type galaxies and the overall merger fraction reasonably 
well, and can therefore be used to obtain a plausible prediction 
of the frequency of dry mergers in a hierarchical universe. The 
predicted dry merger fraction is consistent with the observa- 
tions to within their combined uncertainties. Thus, both obser- 
vations and theory lend strong support to the notion that major 
spheroid mergers are an important driver of the evolution of lu- 
minous early-type galaxies in recent epochs. 
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APPENDIX 

A. ARE FLY-BYS IMPORTANT? 

Tran et al. (2005) recently showed that all luminous dry 
merger candidates in the outskirts of the z ~ 0.8 MS 1054 galaxy 
cluster with sufficient S/N to measure velocity offsets were in- 
deed bound R < 10 kpc pairs (with velocity differences less than 
200 kms" 1 ). Noting that Tran et al. (2005) and van Dokkum et 
al. (1999) argue that most of these galaxies are in the cluster 
infall region (and therefore are in group environments which 
are about to be incorporated into the cluster proper), their study 
supports the notion that the majority of the dry merger candi- 
dates studied in this paper are bound and will likely merge. 

Yet, one of the few published examples of an E-E interac- 
tion in the local universe (NGC 4782/4783) has been shown 
by Madejsky, Bender & Mollenhoff (1991) to have an unbound 
hyperbolic orbit. Based on simulated images from simulations 
of unbound dry mergers similar to those discussed in the main 
text, we find that morphological disturbances generated by fly- 
by interactions are short-lived, <C 50 Myr for systems of a sim- 
ilar luminosity to the ones we observe here. Nonetheless, if the 
number of fly-by interactions is much larger than the number 
of mergers, fly-bys could significantly contaminate a spheroid 
major merger sample. 

Because galaxies with velocity differences greatly in excess 
of their internal velocity dispersion are extremely unlikely to 
merge, it is worth understanding the likely contribution of high- 
velocity fly-bys to tidally-disturbed spheroid-dominated sys- 
tems in a cosmic-averaged sense. Relative velocity difference 
is a strong function of the velocity dispersion of the local en- 
vironment. We explore the likelihood of producing a given 
tidal effect in different environments by adopting the impulse 
approximation. We wish to work out the number of interac- 
tions as a function of environment, characterized by the ve- 
locity dispersion of the local environment ct. The number of 



interactions Ni„ t oc n^^ SQ r 2 , where n E /so is the number density 
of E/S0 galaxies and r is the pericentric distance. The de- 
pendence of cross-section on velocity can be estimated by as- 
suming that the relative velocity v oc ct. In order to work out 
the cross-section for interaction producing a given displace- 
ment in stars s as a function of v, we note that tidal forces 
F oc r~ 3 . Since interaction timescale t oc 1/v, the displace- 
ment in stars s oc Ft 2 oc r~ 3 v~ 2 . Thus, given a constant tidal 
displacement of stars s, the impact parameter r oc v" 2 / 3 . The 
number of E/S0 galaxies can be estimated assuming a constant 
group stellar M/L; « E/S0 oc /E/soMiaioAhaio- Since M hMo oc ct 3 
and rhaio oc ct, « E /so oc f E /so alone; we assume for simplicity 
/e/so oc ct oc v. Thus, Ni„, oc v 2 / 3 . To estimate the fraction of 
galaxies with this given displacement of stars owing to tidal 
affects, /i nt , one divides by the total number of E/S0 galaxies 
Ne/so /E/so-^haio oc v 4 ; thus, fiat oc v" 10 / 3 . A cosmic-averaged 
fraction must be further weighted by the halo mass function of 
group and cluster-sized halos, $(ct) oc ct~ 3 . Thus, despite the 
dominance of E/S0 galaxies in clusters, the vast majority of 
tidally-disturbed spheroid-dominated galaxies are in the low- 
est velocity-dispersion groups capable of hosting E/S0 galaxies, 
and are therefore likely to merge. 

A complementary analysis following Makino & Hut (1997), 
asking what the average environment of spheroid mergers should 
be (using the scalings above), comes to a similar conclusion: 
/mrg oc 1 /ct 2 , and weighting by the mass function of groups and 
clusters, the cosmic-averaged fraction is again completely dom- 
inated by low-speed E/S0 mergers in low-mass groups. 
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